Magnetic fields in cosmological simulations of disk galaxies

Observationally, magnetic fields reach equipartition with thermal energy and cosmic rays in the interstellar medium of disk galaxies such as the Milky Way. However, thus far cosmological simulations of the formation and evolution of galaxies have usually neglected magnetic fields. We employ the moving-mesh code \textsc{Arepo} to follow for the first time the formation and evolution of a Milky Way-like disk galaxy in its full cosmological context while taking into account magnetic fields. We find that a prescribed tiny magnetic seed field grows exponentially by a small-scale dynamo until it saturates around $z=4$ with a magnetic energy of about $10\%$ of the kinetic energy in the center of the galaxy's main progenitor halo. By $z=2$, a well-defined gaseous disk forms in which the magnetic field is further amplified by differential rotation, until it saturates at an average field strength of \sim 6 \mug in the disk plane. In this phase, the magnetic field is transformed from a chaotic small-scale field to an ordered large-scale field coherent on scales comparable to the disk radius. The final magnetic field strength, its radial profile and the stellar structure of the disk compare well with observational data. A minor merger temporarily increases the magnetic field strength by about a factor of two, before it quickly decays back to its saturation value. Our results are highly insensitive to the initial seed field strength and suggest that the large-scale magnetic field in spiral galaxies can be explained as a result of the cosmic structure formation process.Comment: 5 pages, 4 figures, accepted to ApJ

Constrained Simulations of the Magnetic Field in the Local Supercluster and the Propagation of UHECR

Magnetic fields (MF) in the Local Supercluster (LSC) of galaxies may have profound consequences for the propagation of Ultra High Energy Cosmic Rays (UHECR). Faraday rotations measurements provide some informations about MF in compact clusters. However, very few is known about less dense regions and about the global structure of MF in the LSC. In order to get a better knowledge of these fields we are performing constrained magnetohydrodynamical simulations of the LSC magnetic field. We will present the results of our simulation and discuss their implications for the angular distribution of expected UHECR deflections.Comment: 4 pages + 1 figure. Published on the Proceedings of the 28th International Cosmic Ray Conference, Tsukuba, Japan (2003

Helium-ignited violent mergers as a unified model for normal and rapidly declining Type Ia Supernovae

The progenitors of Type Ia Supernovae (SNe Ia) are still unknown, despite significant progress during the last years in theory and observations. Violent mergers of two carbon--oxygen (CO) white dwarfs (WDs) are one candidate suggested to be responsible for at least a significant fraction of normal SNe Ia. Here, we simulate the merger of two CO WDs using a moving-mesh code that allows for the inclusion of thin helium (He) shells (0.01\,\msun) on top of the WDs, at an unprecedented numerical resolution. The accretion of He onto the primary WD leads to the formation of a detonation in its He shell. This detonation propagates around the CO WD and sends a converging shock wave into its core, known to robustly trigger a second detonation, as in the well-known double-detonation scenario for He-accreting CO WDs. However, in contrast to that scenario where a massive He shell is required to form a detonation through thermal instability, here the He detonation is ignited dynamically. Accordingly the required He-shell mass is significantly smaller, and hence its burning products are unlikely to affect the optical display of the explosion. We show that this scenario, which works for CO primary WDs with CO- as well as He-WD companions, has the potential to explain the different brightness distributions, delay times and relative rates of normal and fast declining SNe Ia. Finally, we discuss extensions to our unified merger model needed to obtain a comprehensive picture of the full observed diversity of SNe Ia.Comment: accepted for publication by ApJL, significant changes to first version, including addition of merger simulatio

Stellar feedback by radiation pressure and photoionization

The relative impact of radiation pressure and photoionization feedback from young stars on surrounding gas is studied with hydrodynamic radiative transfer (RT) simulations. The calculations focus on the single-scattering (direct radiation pressure) and optically thick regime, and adopt a moment-based RT-method implemented in the moving-mesh code AREPO. The source luminosity, gas density profile and initial temperature are varied. At typical temperatures and densities of molecular clouds, radiation pressure drives velocities of order ~20 km/s over 1-5 Myr; enough to unbind the smaller clouds. However, these estimates ignore the effects of photoionization that naturally occur concurrently. When radiation pressure and photoionization act together, the latter is substantially more efficient, inducing velocities comparable to the sound speed of the hot ionized medium (10-15 km/s) on timescales far shorter than required for accumulating similar momentum with radiation pressure. This mismatch allows photoionization to dominate the feedback as the heating and expansion of gas lowers the central densities, further diminishing the impact of radiation pressure. Our results indicate that a proper treatment of the impact of young stars on the interstellar medium needs to primarily account for their ionization power whereas direct radiation pressure appears to be a secondary effect. This conclusion may change if extreme boosts of the radiation pressure by photon trapping are assumed.Comment: 18 pages, 19 figures (main results presented in 13 pages, 10 figures; extended appendix for RT tests with extra 9 figures). Accepted for publication in MNRAS after tiny change

Galactic Centre stellar winds and Sgr A* accretion

(ABRIDGED) We present in detail our new 3D numerical models for the accretion of stellar winds on to Sgr A*. In our most sophisticated models, we put stars on realistic orbits around Sgr A*, include `slow' winds (300 km/s), and account for radiative cooling. We first model only one phase `fast' stellar winds (1000 km/s). For wind sources fixed in space, the accretion rate is Mdot ~ 1e-5 Msun/yr, fluctuates by < 10%, and is in a good agreement with previous models. In contrast, Mdot decreases by an order of magnitude for stars following circular orbits, and fluctuates by ~ 50%. Then we allow a fraction of stars to produce slow winds. Much of these winds cool radiatively, forming cold clumps immersed into the X-ray emitting gas. We test two orbital configurations for the stars in this scenario, an isotropic distribution and two rotating discs with perpendicular orientation. The morphology of cold gas is quite sensitive to the orbits. In both cases, however, most of the accreted gas is hot, with an almost constant Mdot ~ 3e-6 Msun/yr, consistent with Chandra observations. The cold gas accretes in intermittent, short but powerful episodes which may give rise to large amplitude variability in the luminosity of Sgr A* on time scales of 10s to 100s of years. The circularisation radii for the flows are ~ 1e3 and 1e4 Rsch, for the one and two-phase wind simulations, respectively, never forming the quasi-spherical accretion flows suggested in some previous work. Our work suggests that, averaged over time scales of 100s to 1000s of years, the radiative and mechanical luminosity of Sgr A* may be substantially higher than it is in its current state. Further improvements of the wind accretion modelling of Sgr A* will rely on improved observational constraints for the wind properties and stellar orbits.Comment: 16 pages, 18 colour figures. Accepted by MNRAS. Full resolution paper and movies available at http://www.mpa-garching.mpg.de/~jcuadra/Winds/ . (v2: minor changes

On the width of cold fronts in clusters of galaxies due to conduction

We consider the impact of thermal conduction in clusters of galaxies on the (unmagnetized) interface between a cold gaseous cloud and a hotter gas flowing over the cloud (the so-called cold front). We argue that near the stagnation point of the flow conduction creates a spatially extended layer of constant thickness $\Delta$, where $\Delta$ is of order $\sim\sqrt{kR/U}$, and $R$ is the curvature radius of the cloud, $U$ is the velocity of the flow at infinity, and $k$ is the conductivity of the gas. For typical parameters of the observed fronts, one finds $\Delta \ll R$. The formation time of such a layer is $\sim R/U$. Once the layer is formed, its thickness only slowly varies with time and the quasi-steady layer may persist for many characteristic time scales. Based on these simple arguments one can use the observed width of the cold fronts in galaxy clusters to constrain the effective thermal conductivity of the intra-cluster medium.Comment: Accepted for MNRAS. 9 pages; 6 b&w figures; 2 colour figure

Back-in-time dynamics of the cluster IE 0657-56 (the Bullet System)

We present a simplified dynamical model of the ``Bullet'' system of two colliding clusters. The model constrains the masses of the system by requiring that the orbits of the main and sub components satisfy the cosmological initial conditions of vanishing physical separation a Hubble time ago. This is also known as the timing argument. The model considers a system embedded in an over-dense region. We argue that a relative speed of $4500 \rm km/s$ between the two components is consistent with cosmological conditions if the system is of a total mass of $2.8\times 10^{15}h^{-1} M_\odot$ is embedded in a region of a (mild) over-density of 10 times the cosmological background density. Combining this with the lensing measurements of the projected mass, the model yields a ratio of 3:1 for the mass of the main relative to that of the subcomponent. The effect of the background weakens as the relative speed between the two components is decreased. For relative speeds lower than $\sim 3700\rm km/s$, the timing argument yields masses which are too low to be consistent with lensing.Comment: 5 pages, 3 figures, submitted to MNRA